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We prove that the two fluid Landau hydrodynamic equations, when applied to a gas interacting 
with infinite scattering length (unitary gas) in the presence ol harmonic trapping, admit exact 
scaling solutions of mixed compressional and surface nature. These solutions are characterized by 
a linear dependence of the velocity field on the spatial coordinates and a temperature independent 
frequency which is calculated in terms of the parameters of the trap. Our results are derived in the 
regime of small amplitude oscillations and hold both below and above the superfluid phase transition. 
They apply to isotropic as well as to deformed configurations, thereby providing a generalization 
of Castin's theorem (Y. Castin, C. R. Phys. 5, 407 (2004)) holding for isotropic trapping. Our 
predictions agree with the experimental findings in resonantly interacting atomic Fermi gases. The 
breathing scaling solution, in the presence of isotropic trapping, is also used to prove the vanishing 
of two bulk viscosity coefficients in the superfluid phase. 
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The thermodynamic behavior of Fermi gases interact- 
ing with infinite scattering length (hereafter called uni- 
tary Fermi gases) are known to exhibit peculiar univer- 
sal properties. This follows from the fact that in the 
unitary limit of uniform configurations the only avail- 
able lengths are the interparticle distance and the ther- 
mal wavelength. Analogously the remaining energies are 
the thermal energy fcsT and the Fermi energy Ep = 
|^(37r^n)^/'^, where fcs is the Boltzmann constant and n 
is the particle density. This peculiarity makes these sys- 
tems extremely interesting because their thermodynamic 
behavior takes a universal character, suitable to explain 
several physical features of quite different systems, like 
resonantly interacting atomic gases and neutron matter. 
Systematic experimental [H-Il] and theoretical (see for ex- 
ample 0-01 and refs therein) efforts have been made to 
explore this universal thermodynamic behavior in the su- 
perfluid as well as in the normal phase. The concept of 
unitarity is not restricted to the Fermi gas and applies to 
Bose gases as well [8l-[l0|. although three-body recombi- 
nations cause a fast instability of these systems, making 
their experimental investigation problematic especially at 
low temperature. 

The unitary Fermi gas is known to exhibit unique fea- 
tures also from the dynamic point of view. In particular 
Castin [ll| has shown that exact scaling solutions are 
available if the system is trapped by an isotropic three- 
dimensional harmonic potential. A remarkable example 
is the occurrence of an undamped radial breathing mode 
oscillating at the frequency 2iOho, where ujho is the oscil- 
lator frequency of the harmonic potential. This result is 
remarkable because it concerns a strongly interacting sys- 
tem and its validity is not restricted to small amplitude 
oscillations. Furthermore it holds exactly at all temper- 
atures irrespective of the collisional regime and the value 
of the mean free path. It can be regarded as the strongly 
interacting and quantum analog of the most famous clas- 



sical result derived by Boltzmann for an ideal gas trapped 
by an isotropic harmonic potential [T2j . 

Similarly to the case of the Boltzmann gas, the uni- 
versality of the scaling oscillation breaks down in the 
presence of a deformed harmonic trap. In this case no 
general exact result is available, unless one considers the 
collisional hydrodynamic regime characterized by the oc- 
currence of fast collisions. This is the case considered 
in the present paper where we prove the existence of a 
class of scaling solutions describing small amplitude os- 
cillations of the gas around equilibrium, characterized by 
a temperature independent value of the collective fre- 
quency. These solutions have been already identified in 
the literature at zero temperature (see @ and references 
therein). Their existence is proven here at all temper- 
atures, both below and above the critical temperature 
for superfluidity. For simplicity we consider an axially 
symmetric trapping potential 
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where uj±_ and uJz are, respectively, the radial frequencies 
while = xi+yj, but our proof can be easily generalized 
to the more general case of tri-axial trapping. 

We start from the two fluid hydrodynamic equations 
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first derived by Landau [T^l , which describe the dynamic 
behavior of a gas in its superfiuid phase. Equation ([2]) 
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is the equation of continuity where n is the density of 
the fluid and j = m(n„v„ + ngV^) is the current density 
where we have introduced the superfluid {ug) and nor- 
mal densities (n„) and the corresponding velocity fields 
v„ and Vs. Equation ([3]) shows that the transport of 
entropy (s is the entropy density) is carried by the nor- 
mal component of the fluid. Equation Q instead fixes 
the equation for the superfluid velocity, governed by the 
chemical potential (p = fJ^iji, T) is the local value of the 
chemical potential, determined by the equation of state 
of uniform matter). Finally Eq. ([S]) is the Euler equa- 
tion for the current fixed by the gradient of the pressure 
P and of the external potential Vext- Above the super- 
fluid critical temperature the superfluid density vanishes 
and Eqs. ([2113]) and ^ reduce to the usual hydrodynamic 
equations of a normal fluid. 

In the following we prove the existence of exact solu- 
tions of the hydrodynamic equations, corresponding to 
the first sound like ansatz v„ = = v for the velocity 
field of the normal and superfluid components with: 



v = /3(i)r_L + <5(t)z 
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Above Tc, where the superfluid density is absent, the 
same anstaz applies to the velocity held of the fluid. In 
Eq. ([6]) we have considered a velocity flow with symmet- 
ric behavior in the x and y directions, corresponding to 
excitations carrying angular momentum along the z-th 
direction = 0. The generalization of the formalism to 
scaling excitations carrying angular momentum £^ = ±1 
and iz = ±2 is straightforward. 

The choice ([6]) is accompanied by the following scaling 
transformations for the density 

„2a(t)+7(t)„p(r') 
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and for the entropy density 



s(r,t) 
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where r' = {e"^ 'x,e°'^ 'y,e'^^ 'z) is the scaled spatial 
variable and no and so are the particle density and en- 
tropy density calculated at equilibrium. The prcfactor 
g2Q(t)-i-7(t) ^jjg above equations ensures the normaliza- 
tion condition for the density and the conservation of 
total entropy. The local entropy per particle s(r)/n(r) 
is also conserved. At unitary the entropy density can 
be written in the form s(n,T) — nfe{T/Tp{n)) where 
Tp = Ep/ks is the Fermi temperature and fe is a uni- 
versal function that can be derived by the knowledge of 
the equation of state (see for example [ij])j but whose ex- 
plicit form is irrelevant for the proof of our theorem. The 
ansatz ([8|) then requires that the ratio T/Tp{n) should 
be conserved by the scaling transformation. This implies 
that the temperature should exhibit the position inde- 
pendent scaling law 



T{t) = e(2"(*)+'r(t))2/32n^ 
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where Tq is the temperature of the gas at equilibrium. 
Finally the chemical potential, which due to dimen- 
sionality arguments can be written in the form /i = 



Tp{n)f^{T/Tp(n)) where is a dimensionless function, 
exhibits the following scaling behavior 



M(r,t)=e(2"(*)+^W)2/Vo(no(r')) 
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with fiQ calculated at the equilibrium temperature Tq. 
The fact that temperature fluctuations associated with 
the scaling solutions are uniform in space represents a pe- 
culiar feature of these collective oscillations. It implies, in 
particular, that as a consequence of the thermodynamic 
relationship VP ~ s VT + nV/i and of the ansatz ^ , the 
equations (|4]) and ([5|) for the superfluid velocity and for 
the current are exactly equivalent. 

We now prove that the above scaling ansatz actually 
corresponds to an exact solution of the hydrodynamic 
equations. From the equation of continuity one finds the 
following relationship 



[(2d + 7) + (2/3 + S)]no 
+ {l3 + a)ri_ ■ V±no 
+ {S + j)zVzno = 



which implies the identities 

a = — /3, 7 = —(5 



(11) 



(12) 



The same conditions permit to satisfy the equation for 
the entropy density. 

Since at equilibrium the chemical potential satisfies the 
condition V/Xo(r) = — mcj^r^ — jtiw^z (see Eq.Q), the 
equation for the superfluid velocity (or, equivalently, the 
equation for the total current) takes the simplifled form 
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By looking for time dependent solutions proportional 
to g-*"^* and expanding the exponentials of Eq. p^ up to 
terms linear in a and 7, one flnally derives the coupled 
equations 



and 
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yielding the temperature independent result 

uj^^(^ + ± i Vl6A4 - 32A2 + 25) ljI (16) 

for the collective frequencies where A = iOz/uj± is the 
aspect ratio of the trap. 

For isotropic trapping (A = 1, i.e. u!± = = uJho) 
the corresponding modes are the uncoupled monopole 
(breathing) mode with uj = 2ojho and the surface 
quadrupole mode with uj = -v/^i^ho- For highly elongated 
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traps (A <C 1, i.e. lu^ ^ uj±) the two solutions are in- 
stead to = a/10/3wj_ and cu = i/l^/Sw^, in agreement 
with the results already derived in the literature at zero 
temperature (see for example [l^ [11] ) . 

The above results provide a generalization of the the- 
orem of where it was shown that for isotropic trap- 
ping the monopolc breathing mode of the unitary Fermi 
gas oscillates with frequency uj = 2oJho, independent of 
temperature. Differently from the case of our results 
hold, however, only in the dissipationless hydrodynamic 
regime and in the limit of small amplitude oscillations. 
The temperature independence of the frequency of the 
scaling solutions of the unitary Fermi gas has been al- 
ready confirmed in experiments (TtI . [18| . 

The scaling solutions discussed above are character- 
ized by temperature variations and by an axial velocity 
field Vz independent of the radial coordinates x and y. 
These are the conditions required, in general, to apply 
the ID hydrodynamic equations in highly elongated con- 
figurations (see [3, 13) which have been recently suc- 
cessfully applied to describe the experimental results of 
p^ . In general these ID like conditions are ensured by 
the effective role of the thermal conductivity and of the 
viscosity which cause the absence of gradients in the ra- 
dial direction and are favored by the presence of a tight 
radial confinement. In the case of the low frequency os- 
cillations considered in the present work the absence of 
radial gradients is automatically ensured by the form of 
the scaling transformation. This explains, in particular, 
why the frequency u) = yJVlJZuJz of the axial breathing 
mode, here derived in a 3D framework in the highly elon- 
gated limit ujz a;_L, coincides with the predictions of 
the ID hydrodynamic equations [lij . 

Let us also mention that the proof of the temperature 
independence of the frequency of the scaling solutions 
can be derived also for the = ±1 and £z = ±2 ex- 
citations characterized by velocity fields proportional to 
V[2;(a; ± iy)] and V(a; ± iy)^, respectively and yielding 



the values i/a 



Jz and y/2uj± for the corresponding 



frequencies. Differently from the iz = solutions dis- 
cussed above the results for the = ±1 and = ±2 
arc not restricted to the unitary Fermi gas, but simply 
require the applicability of the hydrodynamic equations, 
being solutions characterized by divergency free velocity 
fields (surface excitations). They hold in particular for 



both Fermi and Bose gases in the presence of harmonic 
trapping. 

Let us finally discuss a non trivial implication of the 
scaling solutions concerning the behavior of the bulk vis- 
cosity coefficients. According to Khalatnikov [2^ the 
entropy production per unit volume associated with a 
hydrodynamic flow is R/T, where R is the so called dis- 
sipative function in the superfluid phase, defined by 



R = 



( dkv„ 



diVnk 



+2CiV • v„V • mns (v^ - v„) 4- 

+C2 (V • v„)' + Cs [V • mn, (v, - v„)]' (17) 

+ {K/T){WTf . 

In the above equation ij is the shear viscosity, Ci, 
and Cs are the three bulk viscosity coefficients appearing 
in the superfluid phase, while n is the thermal conduc- 
tivity. For the scaling modes discussed in the present 
paper VT ~ and the velocity fields for the normal 
and superfluid components coincide. Furthermore for 
the monopole breathing mode in an isotropic trap one 
has v„ = = /3r and the first term, proportional to 
the shear viscosity, identically vanishes. It follows that 
in this case only the term with (2 survives. However, ac- 
cording to Castin theorem [ll| the dissipation associated 
with the breathing oscillation must be zero and we then 
conclude that C2 = 0. On the other hand the positive- 
ness of R implies that C,f < C2C3 • Thus also Ci must 
vanish, the only surviving bulk viscosity coefficient being 
(^3. Above Tc, where one can introduce only one bulk 
viscosity term, the viscosity coefficient ( should be also 
zero in order to ensure the absence of dissipation. The 
same results, concerning the value of the bulk viscosity 
coefficients, were previously derived by Son [2l| using dif- 
ferent considerations. Actually, our derivation provides 
a simple foundation to Son's heuristic argumentation. 
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